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We study a generalization of the Skyrme model with the inclusion of a sixth-order term and a 
generalized mass term. We first analyze the model in a regime where the nonlinear a and Skyrme 
terms are switched to zero which leads to well-behaved analytical BPS-type solutions. Adding 
contributions from the rotational energy, we reproduce the mass of the most abundant isotopes 
to rather good accuracy. These BPS-type solutions are then used to compute the contributions 
from the nonlinear sigma and Skyrme terms when these are switched on. We then adjust the four 
parameters of the model using two different procedures and find that the additional terms only 
represent small perturbations to the system. We finally calculate the binding energy per nucleon 
and compare our results with the experimental values. 

PACS numbers: 12.39.Dc, ll.10.Lm 

I. INTRODUCTION 

The Skyrme model [lj is nowadays one of the strongest candidates for a description of the low-energy regime of 

QCD. Developed in the beginning of the 60's by T.H.R. Skyrme, it consists of a nonlinear theory of mesons where its 

main feature is the presence of topological solitons as solutions. Each of these solutions is associated with a conserved 

II ■ topological charge, the winding number, which Skyrme interpreted as the baryon number, thus leading him to state 

Oh' that the solitons are baryons emerging from a meson field. The -^- expansion of QCD introduced by t'Hooft |2| in 

:*5 ■ the mid-70's and the later connection from Witten Q with the model developed by Skyrme brought some support to 

\* ' this interpretation. 

Since its original formulation, the Skyrme model has been able to predict the properties of the nucleon within 
a precision of 30%. Several modifications to the model have been considered to improve these predictions, from 
the generalization of the mass term [jjfl to the explicit addition of vector mesons [3, 3, aside form higher order 
terms in derivatives of the pion field [4|. Unfortunately, the analysis of these models has been hampered by their 

^■s' nonlinear nature and the absence of analytical solutions. Indeed, all the solutions rely on numerical computation at 

,—{ • some point whether one uses the rational map ansatz [9|, which turns out to be a rath er g ood approximation of the 
angular dependance, or a full fledge numerical algorithm like simulated annealing (lfj, lll| to find an exact solution 
of the energy functional. Clearly, even a prototype model with analytical solutions would allow to go deeper in the 

p^ investigation of the properties and perhaps identify novel features of the Skyrmions. 

7—1 , In a recent study, Adam, Sanchez-Guillen and Wereszczynski (ASW) [12j obtained an analytical solution by consid- 

ering a model consisting only of a term of order six in derivatives and a potential which correspond to the customary 
mass term for pions in the Skyrme model [13| . Their calculations lead to a compacton-type solution with size growing 
as n 3 , with n is the winding number, a result in general agreement with experimental observations. Another important 
remark on their study is that their solutions are of BPS-type, i.e. they saturate a Bogomolny's bound. Even though 
physical nuclei do not saturate such a bound, the small value of the binding energy may be one of the motivation 
for solution of this type. Let us also mention that recently Sutcliffe [1J| found that BPS-type Skyrmions also emerge 
from models when a large number of vector mesons are added to the Skyrme model. However, the analysis of ASW 
neglects rotational or isorotational energies of nuclei and perhaps the oddest feature of the model is that it does not 
contain any of the terms that Skyrme originally introduced in his model, the nonlinear sigma and so-called Skyrme 
terms which are of order 2 and 4 in derivatives respectively. Being a effective theory of QCD, there is nevertheless 
no reason to omit such contributions. In their work ASW further suggest that their analytical solutions found could 
be used to compute the contributions from the terms of the original Skyrme lagrangian assuming they are small and 
do not affect significantly the overall solutions. Unfortunately the nature of the solution leads to singularities in the 
computation of the energies related to the nonlinear sigma and Skyrme terms. 

In this work, we find analytical BPS-type solutions for a Lagrangian similar to the one in [12j which allows to consider 
contributions from the original Skyrme Lagrangian as small perturbations. The analysis also includes contributions 
for (iso)rotational energies providing a more realistic description of nuclei. The paper is divided as follows: In section 
HI1 we introduce the general form of this generalized Skyrme model and find expressions for the static energies. Next, 
we quantify semiclassically the zero modes of the Skyrmions, which will allow to compute rotational contributions 
to the total energy coming from the spin as well as the isospin of the nuclei. In section IIV| we choose an adequate 
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potential (or mass term) and switch off the nonlinear sigma and Skyrme terms. We then find a simple analytical 
form of the BPS-type solutions for the remaining Lagrangian. It turns out that all the properties of the nuclei can 
be calculated analytically. In section |V] we use the solution to compute the properties of the full Lagrangian. Fitting 
the different parameters of the model with nuclear mass data [15j we verify that the contributions from the nonlinear 
sigma and Skyrme terms remain small and that the analytical solution is a good approximation. 

II. LAGRANGIAN OF THE SKYRME MODEL 

The model proposed by ASW is based on the Lagrangian density 

C = Ce-u?V = -f^Tt ([£„,£„] [L",L A ] [L x ,L,]) - ^ 2 Tr [1 - U] . (1) 

where U — <fio + iTi<fii is the SU(2) matrix representing the meson fields and L M = U^d^U is the left-handed current. 
The model leads to BPS-type solitons. The constants A and /j, are the only free parameters of the model with units 
MeV -1 and MeV 2 respectively. Using scaling arguments, one can show that the term of order 6 in field derivatives, 
Cq, prevents the soliton from shrinking to zero size while the second term, often called the mass term, stabilize the 
solution against arbitrary expansion. 

On the other hand, the original Skyrme model consists of the two completely different terms 

£ = £>nlo + £-Sk (2) 

with 

C nLa = aTr [L^] , C sk = /3Tr ([L M , L„f) (3) 

the nonlinear sigma and so-called Skyrme terms which are of order 2 and 4 in derivatives respectively. Here [a] = MeV 
and (3 is a dimensionless constant. 

We shall consider here a model containing the four terms i.e. an extension of the Skyrme model with a sixth order 
term in derivatives and generalized mass term. The Lagrangian density reads 

C = -fi 2 V(U)-aTT[L^}+f3Tv{[L^L u } 2 ) - | A_Tr ([L M , L v ] [L V ,L X ] [L X ,L„]) (4) 

We are interested in the regime where a and (3 are small so that C n i a and Csk can be considered as small perturbations 
to ((TJ. Usually, the potential V(U) is chosen such that it reproduces the mass term for pions when small fluctuations 
of the fields are considered 

U = e»r*%„ 1 + 2^^-. ( 5 ) 

where F^ = 4y / a is the pion decay constant. Since U is an SU (2) matrix, the meson fields obey the condition 

0o + 0? = 1 (6) 

to limit the degrees of freedom to three. The boundary condition 

[/(r^co)=I 2x2 (7) 

with X 2X 2 the two-dimensional unit matrix, ensures that each solution for the Skyrme field falls into a topological 
sector characterized by a conserved topological charge 

B = ~~S^ j " dVrr (M^^D- (8) 

The static energy can then be calculated using 

E sta t = - fd 3 x£ st at- (9) 



We may conveniently write a general solution as 



where n is the unit vector 



or 



U = e m ' Tf = cos F + in • r sin F 



n = (sin cos $, sin 9 sin $, cos 0) 



4> a = (cos F, sin F sin cos $, sin F sin sin <J>, sin F cos 0) . 
Following ASW [13, we consider solutions of the form that saturates the Bogomolny's bound for (JTJ) 

F = F(r), = 0, $ = n0 

where n is an integer. The static energy @ becomes 
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where F' = dF/dr and the topological charge is simply B = n. 

In order to represent physical nuclei, we have to quantize the solitons using a semiclassical method described in the 
next section. Using the appropriate spin and isospin numbers, we will then be able to calculate the total energy for 
each nuclei. 



III. QUANTIZATION 



Because the topological solitons occupy a spatial volume that is nonzero, usual quantization procedures are no 
longer available. We therefore have to use a semiclassical quantization method by adding an explicit time dependence 
to the zero modes of the Skyrmion. Performing time-dependent (iso) rotations on the Skyrme field by SU(2) matrix 
A(t) and Bit) yield 



U (r, t) = A{t)U (R(B(t))r) A(t) 



(13) 



where Rij(B(t)) = ^Tr [tjStjjB'] is the associated 5*0(3) rotation matrix. Upon insertion of this ansatz in the 
time-dependent part of Q, we write the rotational lagrangian as 



C-rot = -^aiUijCLj - QiWijbj + -biVijbj 



with Uij, Vij and Wij the inertia tensors 



Uij = - f d 3 a; J2aTr (T % T 3 ) + 4/Hr ([L k ,T$ [L k ,T 3 ]) 

9A 2 1 

+—Tr([T l ,L k ] [L k ,L n ] [L n ,Tj]) j, 



(14) 



(15) 



Vij = -e ik ie jmn / d xx k x m \ 2aTr (LiL n ) + 4/3Tr (\L P ,L{\ [L p ,L n ]) 
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Wij = e jk i 
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f d^xx,, UaTr (T^) + 4/3Tr ([L p , Tj] [L p , L r 



2 Tr ([Tj, L m ] [L m , L n ] [L„, L{\) 



(17) 



and Ti = iW [^,[/]. Assuming a solution of the form ([TU]) . the inertia tensors becomes all diagonal and furthermore, 
one can show that (7 n = C/ 2 2 7^ E^33 with similar identities for the Vij and Wij tensors. Finally the general expressions 
for the moments of inertia coming from each pieces of the Lagrangian read 
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(18) 
(19) 



and the general expression for t/33 can be obtained by setting n = 1 in the integrand of (|18|) . It turns out that 
expressions (fT5)) - (fT7)) leads to W\i = W22 = for \n\ > 2 and n 2 J7 33 = n^ = V33. Otherwise, for \n\ = lwe have 
spherical symmetry and, 
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Following Houghton and Magee [Tl|, we now write the rotational hamiltonian as 
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with (Ki) Li the body-hxed (iso)rotation momentum canonically conjugate to a% and hi respectively. The expression 
for the rotational energy of the nucleon has been obtained in [ll| and reads, for a spherical symmetry 
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For the deuteron, the rotational energy has been calculated assuming an axial symmetric solution |16| 
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E rot = 
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which reduces to 
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for the axial ansatz (fTU|) . It is easy to calculate the rotational energies for nuclei with winding number \n\ > 3. The 
axial symmetry of the solution imposes the constraint L 3 + nK 3 = which is simply the statement that a spatial 
rotation by an angle 9 about the axis of symmetry can compensated by an isorotation of —n6 about the T3 axis. It also 
implies that n 2 U 33 = nW 33 = V33. Recalling that W\\ = W21 = for these values of n, the rotational hamiltonian 
reduces to 
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These momenta are related to the usual space-fixed isospin (I) and spin (J) by the orthogonal transformations 

It = -RiAJijKj, (26) 



A = -R(A 2 )J 3 L r 

According to (|26|) and (|27|) , we see that the Casimir invariants satisfy K 2 
the rotational hamiltonian is given by 
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IV. BPS-TYPE SOLUTIONS 

Let's consider a model similar to |12| composed of the term of order six in derivatives plus a potential by setting 
a, /3 = 



3 A 



2 



C = -oT^^{\L^L v ][L\L x ] [L X ,L^])-^V(U). (29) 
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Using the results of section UH the static energy is 

f /9A 2 n 2 sin 4 F ,, , , \ 
E stat =i7r I dr(— ^-F' 2 + ix 2 V(U)J . (30) 

The minimization of the static energy of the soliton, leads to the differential equation for F 

9A2 ^ Sin %f^W^ = 0. (31) 



4 2r 2 
A change of variable z — ' 9 2f f" allows to rewrite (|3L[) in a simple form 



This last equation can be integrated 
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sin 2 Fd z [sin 2 F (d z F)] - M 2 — = 0. (32) 



i sin 4 F(F Z ) 2 = V (33) 



and inserting the expression for z, provide an expression which amounts to a statement of equipartition of the energy, 
i.e. the term of order 6 in derivatives and the potential contribute equally to the total energy. ASW has shown that 
a solution of (|3"3"|) saturates the Bogomolny's bound [12|- From (|3"3"|) we obtain the following useful relation between 
the function F and the potential 

sin 2 F 
dF^= =±(z- z ) (34) 

with zq an integration constant. 

Now comes the time to choose a specific potential. The choice for mass term of the Skyrme is not unique and 
indeed has been the object of several discussions J4|, la LL3| • The usual mass term V = 1 — cos F was considered in 
G3- Solving dMJ) for F 

„, . I 2arccos (vr 3 ) for re 0,^~3 ,„_. 

F(r) = i v i L_ A J (35) 

^0 for r > v 3 

where v = 18 M , is a constant depending on the parameters A, \x and n. Note that F' diverges as r — > v~^. Since this 
solution saturates the Bogomolny's bound, the static energy is proportional to the baryon number B = n . 

A question arises as how would the nonlinear a and Skyrme term affect the energy of such Skyrmions. Switching 
them on slowly by moving a and (3 away from zero could give an estimate of their contributions. Unfortunately, it 
turns out that simply substituting the solution (|35[) in the expression for energy associated with the full Lagrangian 
leads to divergences. So, however small the parameters a and (3 are, these BPS solutions cannot be considered as 
appropriate approximations of the solutions for ([J| . 

Yet it could be interesting to analyze the full Lagrangian (0| in a regime close to a BPS Skyrmion. For this purpose, 
we propose to write the potential in the form of the generalized mass term introduced by [4j 

oo 

V = - Y^ CfcTr [U k + U^ - 2] 
fc=i 

= -4 7 r£y'r 2 dr8C fc sin 2 ^. (36 ) 
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The main motivation for this choice is that the potential can be written in a simple form in terms of pion fields. 
Furthermore this particular framework insures that one recovers the chiral symmetry breaking pion mass term — ^m^ir- 
it in the limit of small pion field fluctuations provided 

E fc2 ^ = - ^ (37) 

fc=l 

For practical purposes, one requires (i) an expression of the potential that is simple enough to allow the analytical 
integration of the left-hand side of equation (|34p . (ii) that the results leads to an invertible function to be able to 
write the chiral profile F as a function of r and finally (iii) that F (r) is well behaved. A most convenient choice is 

^rf^jcos 6 ^), (38) 

Expanding the expression (|55)) , the coefficients C& are 

LI Li Li LI 

Cl = ~128' C2= 128' a3= 128' C4= 512' ° k>A = °' (39) 

Integrating (|3"4")) we get the general solution 

F(r) = 2 arccos U^^-^A (40) 

with v = 18 M . . In order that the baryon number corresponds to \B\ =n, one must require that 

F(oo) - F(0) = Ttt 

for B positive or negative respectively. Accordingly we choose the boundary conditions F(0) = and F(po) = T 71 " 
which sets the integration constant r$ to zero and allows to write 

F{r) = t2 arccos (e" 1 "' 3 ) (41) 

where we use the absolute value to dispose of the sign ambiguity of the arccos function. Note that here contrary 
to [12j we do not get a compacton type solution but a well behaved function with a continuous first derivative. All 
calculations regarding energy can be performed analytically, i.e. static energy and the moments. For example, the 
baryon density is given by the radial function 
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which upon integration leads to baryon number B = n. Experimentally the size of the nucleus is known to behave as 

R = R B^ = (1.25 fin) £3 

where Rq = 1.25 fin. It is interesting to note that the baryon number distribution is zero at r = but has maximum 
value g -/^ independent of n which is positioned at 

r m ,x=(—bx(J\Y Bi (42) 

where r max here is in units of MeV -1 . Accordingly the size of the nucleus r max is proportional to Bs with Rq 
depending only on the ratio X/fi. Similarly expressions can be obtained for energy and moment of inertia density. 
Using g0|), they yield 

E = 2mr[i\ 

«—"^-^Tiy l '»- J '(S) , ' ft ' (!)("■ »'-••»') ' 43 » 

where T (x) is the gamma function. Combining these results in 



H r 



2U n 



J(J+1) 



3n 2 



An 2 



1(1 + 1) 



Kl 



(44) 



Note that this last result only hold for a = /3 = and the solution (|40[) . The last term is either zero or negative. 
Depending on the dimension of the spin and isospin representation, the diagonalization of this hamiltonian will lead 
to a number of eigenstates. We are interested in the lowest eigenvalue of H rot which points towards the eigenstate 
|z, Z3 , /C3) \j,J3,h) with the largest possible eigenvalue k%. Since K 2 = I 2 and L 2 = J 2 , the state with highest weight 
is characterized by k 3 = I and 1$ = j and since since nuclei are build of B fermions j < B/2. On the other hand the 
axial symmetry of the solutions implies that k 3 = —l 3 /n. We recall that these solutions are approximations. Then for 
even B nuclei, the integer part of |?3/n.| is 



< Ifc, 



so it leads to |fc3J = 0. Similarly for half-integer spin nuclei, 



So we shall assume for simplicity 



The rotational energy is given by 
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It remains to fix the values of the parameters A and (i. In order to do so, we choose as input parameters the 
experimental mass of the nucleon and for simplicity, a nucleus X with zero (iso)rotational energy (i.e. a nucleus with 
zero spin and isospin) . The total energy of these two states are 
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Solving for A and \i we get 



Ex = E\ n=B = 2B7T/XA 
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As an example, we choose the nucleus X to be the helium-4, the first doubly magic number nucleus. The mass of the 
helium-4 nucleus has no (iso)rotational parts since it has zero spin and isospin. Setting the mass of the nucleon as the 
average mass of the proton and neutron i.e. En = 938.919 MeV and the mass of the helium nucleus to En e — 3727.38 
MeV, we obtain the numerical value A = 0.00491505 MeW 1 and /x = 30174.2 MeV 2 . We shall refer to this set of 
parameters as Set la. 

Experimentally the size of the nucleus is known to behave as 



R = RoB* 



with Rq = 1.25 fm We get a similar behaviour for r max hi (|42p, 

r max = (1.4798 fm) Si (49) 

Combining (|4"B")) with (fT2")l and (|¥5j) , the mass of any nucleus can be expressed as a analytical function of the input 
parameters E^ and Ejj e . In general it depends on the baryon number as well as the spin and the isospin of the 
isotope. The spin of the most abondant isotopes are known. The isospins are not so well known so we resort to the 
usual assumption that the most abundant isotopes correspond to states with lowest isorotational energy, i.e. states 
where the isospin I has the lowest value that the conservation of the third component of isospin I3 allows. Accordingly, 
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Table I shows the results for the a few isotopes. The resulting predictions are accurate to 0.3% or better even for 
heavier nuclei which is rather surprising since the model involves only to two free parameters A and fj.. 

The computation were repeated using as input parameter X = 16 and 40 Ca, two other doubly magic nuclei 
(also shown in Table I, Set lb and Set Ic respectively). These set the parameters to A = 0.00449295 MeV -1 and 
H = 32977.0 MeV 2 and to A = 0.00426504 MeV" 1 and \i = 34717.8 MeV 2 respectively. Using these heavier elements 
as input parameters changes slightly the overall predicting accuracy. Whereas the best overall accuracy is achieved 
using 16 parametrization in Set lb, the lightest isotopes are best described by choosing 4 He as input (Set la). Note 
that the lightest nuclei have lower moments of inertia and get relatively large rotational contribution to their mass. 
Consequently their masses are expected to be more sensitive to the parameters affecting rotational energy. Likewise, 
since the ratio \//j, decreases for X = 16 and 40 Ca, the size of the nucleus also decreases with i? = 1.3943 fm and 
1.3470 fm respectively. 







Table I: Nuclear masses 


(MeV) 
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Nucleus 


Set la 

£ B (N+ 4 He) 
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E B {N+ 16 0) 
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938.919 
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2 H 
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3 H 


2798.13 
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2794.23 
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4 He 
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3723.77 


3721.47 


3727.38 
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6 Li 


5592.02 


5586.73 
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7 Li 


6524.66 


6518.57 


6514.68 


6533.83 
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9 Li 


8387.75 
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8392.75 
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20 Ne 


18636.9 


18618.8 


18607.3 


18617.7 


40 


40 Ca 


37273.8 


37237.7 
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37214.7 


56 


56 Fe 


52183.4 


52132.9 


52100.7 


52089.8 


238 


238TJ 


221780 


221565 


221429 


221696 



Given this unexpected success, one may wonder how switching on the nonlinear a and Skyrme terms can improve or 
affect these results. Indeed the last results suggest that these contributions need not be be very large. This aspect is 
analysed in the next section. 



NONLINEAR a AND SKYRME TERMS 



Let us now consider the full Lagrangian in ^j assuming that the contribution the nonlinear a and Skyrme terms 
can be set arbitrarily small so that (|4"0")) represents a good approximation to the exact solution. Inserting the solution 
in (|12p and in the expression for the various moments of inertia, one get additional contributions proportional to a 
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and as above W\\ = 5 n ,iUn otherwise W\\ = W22 = for |n| > 2 . Again due to the axial symmetry of the ansatz, 



V\\ = U22 ^ U33 while non diagonal elements of Uij are zero. Similar identities also holds for the Vi 



ZJ and Wij tensors. 



Furthermore we have n 2 U33 = 11^33 = V33. Relations (|5H55p bring a clear understanding of the dependence of the 
masses of the nuclei on the various parameters B = n, /x, a, /3 and A as long as a and /3 remain relatively small. 

In order to estimate the magnitude of the parameter a and /3 in a real physical case, we perform two more fits: Set 
II optimizes the four parameters /1, a, (3 and A to reproduce the best fit for the masses of the nuclei and Set III is done 
with respect to the ratio of the binding energy over atomic number, B.E./A. More precisely, we use only a subset of 
table of nuclei [15| composed of the most abondant 144 isotopes (see Fig. 1). This is compared to Set I which was 
determined in the previous section using the masses of the nucleon and 4 He and assuming a = j3 = 0. The results are 
presented in Fig. 1 in the form of B.E./A as a function of the baryon number for Sets la, II, III and experimental 
values. The optimal values of the parameters are 



Table II: Value of parameters for different fits 



Nucleus 



ix (MeV 2 ) 
a (MeV 2 ) 
(dimensionless) 
A (MeV" 1 ) 



Set la 

(N+ 4 He) 



30174.2 





0.00491505 



Set II 

(Masses) 



29841.2 
0.00830341 

-5.48285 x 10~ 7 
0.00496265 



Set III 

(B.E./A) 



29475.7 

0.0316869 

-4.01085 x 10~ 7 

0.00503994 



As suspected the new sets of parameters are very close to Set la. . The nonlinear a and Skyrme parameters a et /3 

1 
are very small but in order to compare, it is best to rescale the static energy with the change of variable u = (4/i/3A) 3 r 



10 



such that the relative weight of each term is more apparent. Then the static energy takes the form 

sin 2 Fl 



*-* = **{¥) J **> ( v + 2 °{^ 



F' 



(n 2 + 1) 



16/3 



16 
9X 2 fi 



sin 2 F 



(n 2 + 1) F' 2 + n 2 



sin 2 F 



9 „,9sin F 
n 2 F' 2 ;— 



(56) 

(57) 



where F' = dF/du and the energy can be expressed in units of -^. For example for Set II (Set III), the nonlinear 
10 
9 




B.E./A 

(MeV) 



o Set la 



o Set III 



a Set II 



' Exper. data 



50 



100 150 

Baryon number 



200 



250 



FIG. 1: Ratio of the binding energy (B.E.) over the atomic number A (or baryon number) as a function of A. The experimental 
data (black squares) are shown along with predicted value for parametization of Set la (empty circles), II (empty triangles) 
and III (empty diamonds) respectively. 



a term is proportional to a I s^~- 



= 3.73524 x 1(T 7 (1.43418 x 1CT 6 ) and the Skyrme term to /3 ( 



16 
9AV 



(—9.2355 x 10 6 ) while the remaining terms are of order one. Furthermore the overall factor I ^H j 



-1.3263 x 10~ 6 

remains approximately the same for all the sets. Looking at the numerical results, we observe nonetheless that these 
two terms are responsible for corrections of the order of 0.01%. Clearly, the small magnitude of these contributions 
provides support to the assumption that (|4"0")) is a good approximation to the exact solution. 
Comparing Set II to the original Skyrme Model with a pion mass term, we may identify 



F n = 4V^ = 0.364474 MeV 

e 2 = ^ 
32/3 

2^ 



-57019 
= 231591 MeV 



(Experiment: F„ = 186 MeV) 
(e = 4.84 for massive pion Skyrme Model) 



/' 



(Experiment: m^ = 138 MeV) 



Set III leads to similar values for F^ , e 2 and to t which are orders of magnitude away for the usual values obtained for 
the Skyrme model. Of course here the nonlinear a and Skyrme terms do not play a significant role in the stabilization 
of the soliton. Indeed the Skyrme term even have the wrong sign so it would destabilize the soliton against shrinking if 
it was not for the contribution of order six in derivatives. The size of the soliton is instead determined by the relative 
magnitude of /x and A so there is no need for F„ and m^ to be close to the nucleoli mass scale as for the original Skyrme 
Model. Perhaps the explanation for such a departure is that the parameters of the model are merely bare parameters 
and they could differ significantly from their renormalized physical values. We note also that the B > 1 solutions of 
the Skyrme Model display a totally different structure compared to the BPS-type solution analyzed here. It is well 
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known that the lowest-energy B = 2,3,4... solutions of the Skyrme Model exhibit respectively toroidal, tetrahedral, 
cubic,... baryon density configurations. Such solution are conveniently represented by an ansatz based on rational 
maps Q. The model at hand here leads to spherically symmetric baryon density at least in the regime of small a 
and j3 where solution (|4"0)) can apply. So it seems that the regime dominated by the /j, and A terms leads to spherical 
configurations whereas the regime dominated by the nonlinear a and Skyrme terms shows totally different baryon 
density distributions. In the absence of a complete analysis, we can only conjecture that the change in configuration 
is related to which of the four terms are responsible for the stabilization of the soliton and at some critical values of 
the parameters there is a transition between configurations. 

Let us now look more closely at the numerical results presented in Fig. 1. These are in the form of the ratio of 
the binding energy (B.E.) over the atomic number A as a function of A which corresponds to the baryon number. 
The experimental data (black squares) are shown along with predicted value for parametization of Set la (empty 
circles), II (empty triangles) and III (empty diamonds) respectively. Clearly Set la is less accurate when it comes to 
reproduce the full set of experimental data but is somewhat successful for the lightest nuclei. This to be expected 
since the fit relies on the masses of the nucleon and 4 He. Yet all predicted nuclear masses are found to be within a 
0.3% precision. In fact the ratio B.E./ A is rather sensitive to small variation of the nuclear masses so the results in 
general are surprisingly accurate. On the other hand Set II, based on the nuclear masses, overestimates the binding 
energies of the lightest nuclei while it reproduces almost exactly the remaining experimental values. The least square 
fit based on B.E. /A, Set III, is the best fit overall but in order to better represent the features of lightest nuclei, it 
abdicates some of the accuracy found in Set II for B > 30. 

This apparent dichotomy between the description of the two regions B < 30 and B > 30 may find an explanation 
in the (iso)rotational contribution to the mass. Indeed light nuclei have smaller sizes and moments of inertia so that 
their rotational energy contributes to a larger fraction of the total mass since the spins and isospins remain relatively 
small. On the other hand the size of heavy nuclei grows as B' and their moments of inertia increase accordingly. 
The spin of the most abundant isotopes are relatively small while isospin can have relatively large values due to 
the growing disequilibrium between the number of proton and the number of neutron in heavy nuclei (see eq. |50"|) . 
Despite these behaviors, our numerical results show that the rotational energy is less that 1 MeV for B > 10 for any 
of the Sets considered and its contribution to B.E./ A decreases rapidly as B increases. On the contrary for B < 10 
the rotational energy is responsible for large part of the binding energy which means that B.E. /A should be very 
sensitive to the way the rotational energy is computed. In our case we approximated the nucleus as a rigid rotator. 
One may argue that if rotational deformations due to centrifugal effects were to be considered, it would lead to larger 
moments of inertia and lower rotational energies. This would predominantly affect the binding energy of the lightest 
nuclei since this is where rotational energy is most significant. Allowing for such deformation would in general require 
the full numerical computation of the solution. An easier way to check for deformation is by allowing the ratio of the 
parameters a = fi/X in the solution (|40[) to vary independently from the /i and A in the model |4|) and by repeating the 
fit with respect to five parameters instead of the four previous ones. This procedure allows for a further adjustment 
of the size of the soliton in terms of a with respect to a given choice of model parameters fi, a, (3 and A and would lead 
to partial deformation of the solution. Such a parametrization is expected to increase both the size and the moments 
of inertia of the soliton and decrease the total mass of the lightest nuclei which would be an improvement over the 
four parameters fit. We evaluated such correction for the nucleon whose relative contribution to mass from rotation 
is the largest using the parameters of Set II and we obtained a modest decrease of the mass of the order of 0.16%. 
Since the rotational energy accounts for much less than 1% of the total energy in most of the nuclei, deformations are 
not generally expected to be very significant. 

VI. CONCLUSION 

We have proposed a 4-terms model as a generalization of the Skyrme Model. In the regime where two of the terms 
are negligible, i.e. a — (3 — 0, we find well-behaved analytical solutions for the static solitons. These saturate the 
Bogomolny's bound with consequence that the static energy is directly proportional to the baryon number B. They 
differ from those obtained by ASW in an important way: their model lead to compactons at the boundary of which 
the gradient of the solution is infinite and so the solution could not be used to approximate the energies in the regime 
where a, (3 ^ 0. Furthermore, one of the major feature of our model is that the form of the solutions allows to compute 
analytically the static and rotational energy and express them as a function of the model parameters and B. Fixing 
the remaining parameters of the model, /i and A, leads to rather accurate predictions for the mass of the nuclei. 

We then used these BPS-type solutions to compute the mass of the nuclei in the regime where a and (3 are small 
but not zero. Indeed fitting the model parameters to provide the best description of the nuclear mass data leads to 
that particular regime where the values of a and /3 turn out to be very small. Yet we find a noticeable improvement 
in the size and B.E. /A predictions with respect to those for the a = f3 = regime. Even though our 4-term model 
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can be considered a simple extension of the massive pion Skyrme Model (different mass term and an additional term 
with six derivatives in pion fields) the solution leads to spherically symmetric baryon densities as opposed to more 
complex configurations for B > 1 standard Skyrmions (e.g. toroidal, tetrahedral. cubic...). These results suggest that 
nuclei could be considered as near BPS-Skyrmions. 

This work was supported by the National Science and Engineering Research Council of Canada. 
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